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Extensive equilibrium Monte Carlo simulations are performed for a three-dimensional Heisenberg 
spin glass with the nearest-neighbor Gaussian coupling to investigate its spin-glass and chiral-glass 
orderings. The occurrence of a finite-temperature chiral-glass transition without the conventional 
spin-glass order is established. Critical exponents characterizing the transition are different from 
those of the standard Ising spin glass. The calculated overlap distribution suggests the appearance 
of a peculiar type of replica-symmetry breaking in the chiral-glass ordered state. 



While experiments have provided convincing evidence 
that spin-glass (SG) magnets exhibit an equilibrium 
phase transition at a finite temperature, the true nature 
of the experimentally observed SG transition and that of 
the low-temperature SG phase still remain open problems 
|jj . A simple Ising model has widely been used as a "re- 
alistic" SG model in the studies, e.g., of the critical prop- 
erties of the SG transition, or of the issue of whether the 
SG state exhibits a spontaneous replica-symmetry break- 
ing (RSB). One should bear in mind, however, that the 
magnetic interactions in many SG materials are nearly 
isotropic, being well described by an isotropic Heisen- 
berg model, in the sense that the magnetic anisotropy is 
considerably weaker than the exchange interaction. In 
apparent contrast to experiments, numerical simulations 
have indicated that the standard spin-glass order occurs 
only at zero temperature in the three-dimensional (3D) 
Heisenberg SG |^-|(|. Although the magnetic anisotropy 
inherent to real materials is often invoked to explain this 
apparent discrepancy with experiments, it still remains 
puzzling that no detectable sign of Heisenberg-to-Ising 
crossover has been observed in experiments which is usu- 
ally expected to occur if the observed finite-temperature 
transition is caused by the weak magnetic anisotropy 

In order to solve this apparent puzzle, a chirality mech- 
anism of experimentally observed spin-glass transitions 
was proposed by one of the authors 0. This scenario 
is based on the assumption that an isotropic 3D Heisen- 
berg SG exhibits a finite-temperature chiral-glass transi- 
tion without the conventional spin-glass order, in which 
only spin-reflection symmetry is broken with preserving 
spin-rotation symmetry. Chirality is an Ising-like multi- 
spin variable representing the sense or the handedness of 
the noncollinear spin structures induced by spin frustra- 
tion. In this scenario, all essential features of many of 
real SG transitions and SG ordered states should be de- 
termined by the properties of the chiral-glass transition 
and of the chiral-glass state in the fully isotropic system, 
while the role of the magnetic anisotropy is secondary 
which "mixes" the spin and the chirality and "reveals" 



the chiral-glass transition as an anomaly in experimen- 
tally accessible quantities. 

Numerical studies on the 3d XY spin glasses have given 
strong support to the occurrence of a finite-temperature 
chiral-glass transition ]^-|l0|]. In the Heisenberg case, 
while previous numerical works agreed in that the stan- 
dard SG order occurred only at T = 0, the question 
whether there really occurs a finite-temperature chiral- 
glass order has remained inconclusive J3||| . Very recently, 
an off-equilibrium Monte Carlo simulation by one of the 
authors has given evidence for the occurrence of a finite- 
temperature chiral-glass order in the 3D Heisenberg SG 
H . However, the full critical properties of the chiral-glass 
transition as well as the properties of the chiral-glass or- 
dered state itself, particularly the question of the possible 
RSB, still remains largely unclear. In this Letter, we per- 
form extensive equilibrium Monte Carlo simulations of 
a 3D Heisenberg SG in order to determine the detailed 
static and dynamic critical properties and to clarify the 
nature of the chiral-glass ordered state. 

Our model is the classical Heisenberg model on a sim- 
ple cubic lattice with N = L 3 spins denned by the Hamil- 
tonian, 
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(1) 



where Si =(Sf , Sf , Sf) is a three-component unit vec- 
tor, and the sum runs over all nearest-neighbor pairs. 
The interactions Jij are random Gaussian variables with 
zero mean and variance J. The local chirality at the i-th 
site and in the fi-th direction, Xn>- is defined for three 
neighboring spins by the scalar 



Xi 



(Si X Si_g fl ), 



(2) 



where e M (p = x, y, z) denotes a unit lattice vector along 
the /x-axis. 

Monte Carlo simulation is performed based on the ex- 
change MC method, sometimes called "parallel temper- 
ing" , which turns out to be an efficient method for ther- 
malizing systems exhibiting slow dynamics [Q. By mak- 
ing use of this method, we have succeeded in equilibrating 
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the system down to the temperature considerably lower 
than those attained in the previous simulations [g|4| . We 
run two independent sequences of systems (replica 1 and 
2) in parallel, and compute an overlap between the chiral 
variables in the two replicas, 



q x ~ Q AT 
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In terms of this chiral overlap q x , the Binder ratio of the 
chirality is calculated by 



9cg(L) 
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where (■ • •) represents the thermal average and [• • •] repre- 
sents the average over bond disorder. For the Heisenberg 
spin, one can introduce an appropriate Binder ratio in 
terms of a tensor overlap q^ v (^-j; = x,y,z) which has 
3 2 = 9 independent components 



1)q(2) 



(fj,,u = x,y,z), 



via the relation, 
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The lattice sizes studied are L = 6, 8, 10, 12 and 16 with 
periodic boundary conditions. Equilibration is checked 
by monitoring the stability of the results against at least 
three-times longer runs for a subset of samples. Sample 
average is taken over 1500 (L = 6), 1200 (L = 8), 640 
(L = 10), 296 (L = 12) and 136 (L = 16) independent 
bond realizations. Note that in the exchange MC simu- 
lations the data at different temperatures are correlated. 
Error bars are estimated from statistical fluctuations over 
the bond realizations. 

The size and temperature dependence of the Binder 
ratios of the spin and of the chirality, gsG and gcG, are 
shown in Fig. 1(a) and (b), respectively. As can be seen 
from Fig. 1(a), gsG constantly decreases with increasing 
L at all temperatures studied, suggesting that the con- 
ventional spin-glass order occurs only at zero tempera- 
ture, consistent with the previous results [^[j6|. Fig. 1(a) 
reveals that gsG for larger lattices (L — 10, 12, 16) ex- 
hibits an anomalous bending around T/J ~ 0.15, sug- 
gesting a change in the ordering behavior in this temper- 
ature range. As can be seen from Fig. 1(b), the curves of 
gcG for different L tend to merge for larger L in the tem- 
perature range where the curves of gsG exhibit an anoma- 
lous bending. Furthermore, on increasing L, the merging 
points gradually move toward higher temperatures, sug- 
gesting that the chiral-glass transition indeed occurs at a 
finite temperature. Since gcG for different L do not cross 
here at gcG > 0, however, it is not necessarily easy to un- 
ambiguously locate the chiral-glass transition point, or 



even to completely rule out the possibility of only a zero- 
temperature transition with rapidly growing {e.g., ex- 
ponentially growing) correlation length. Meanwhile, the 
calculated gcG shows a negative dip whose depth grad- 
ually increases with increasing L, whereas its position 
gradually shifts toward lower temperature. Note that, 
if the systems would not exhibit any finite-temperature 
transition, gco,{T;L — > oo) should be equal to zero at 
any T > and to unity at T = (for the nondegen- 
erate ground state as expected for the present Gaussian 
coupling). Since gcG{T] L — > oo) is a single- valued func- 
tion of T, the existence of a negative dip growing with L 
hardly reconciles with the absence of a finite-temperature 
transition so long as this tendency persists, and suggests 
the chiral-glass transition occurring at T — Tqg > 
at which gec takes a negative value unlike the standard 
cases of the 3D Edwards- Anderson (EA) Ising model or 
the infinite-range Sherrington-Kirkpatrick (SK) model. 
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FIG. 1. Temperature and size dependence of the Binder 
ratios of the spin (a) and of the chirality (b). 



More unambiguous estimate of Tqq can be obtained 
from the equilibrium dynamics of the model. Thus, 
we calculate both the spin and chirality autocorrelation 
functions defined by 



C s (*) = ^E[( S '(*o)-S i (t + to))], (7) 
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C x( f ) = 52[{xMXi»(t + to)}}, (8) 
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where MC simulation is performed according to the stan- 
dard heat-bath updating here. The starting spin con- 
figuration at t = to is taken from the equilibrium spin 
configuration generated in our exchange MC runs. 

Monte Carlo time dependence of the calculated C s (t) 
and C x (t) are shown in Fig. 2 on log- log plots at several 
temperatures for L = 16. We found no significant differ- 
ence in the data of L=12 and 16, and the finite-size effect 
is negligible in our time window. As can be seen from 
Fig. 2(a), C s (t) shows a downward curvature at all tem- 
perature studied, suggesting an exponential-like decay 
characteristic of the disordered phase, consistently with 
the absence of the standard spin-glass order. In sharp 
contrast to this, C x (t) shows either a downward curva- 
ture characteristic of the disordered phase, or an upward 
curvature characteristic of the long-range ordered phase, 
depending on whether the temperature is higher or lower 
than T/J ~ 0.16, while just at T/J ~ 0.16 the linear be- 
havior corresponding to the power-law decay is observed: 
See Fig. 2(b). Hence, our dynamical data indicates that 
the chiral-glass order without the standard spin-glass or- 
der takes place at Tqg/J — 0.160 ± 0.005, below which 
a finite chiral EA order parameter > develops. 
From the slope of the data at T — Tqg, the exponent 
A characterizing the power-law decay of C x (t) ~ t~ x is 
estimated to be A = 0.193 ±0.005. The estimated Xcg is 
in good agreement with the previous estimate of Ref . || , 
T CG /J = 0.157 ±0.01. 



The behavior of the chiral-glass order parameter, or 
the associated chiral-glass susceptibility xcg = 3iV[(g^)], 
turns out to be consistent with this. In the inset of 
Fig. 3, we show the reduced chiral-glass susceptibility 
Xcg = Xcg/x 4 , normalized by the amplitude of the lo- 
cal chirality x 2 = (1/3JV) Y^i,u[(Xifj}]> versus the reduced 
temperature t =\ (T—Tgg)/Tcg I on a log-log plot. From 
the asymptotic slope of the data, the susceptibility ex- 
ponent is estimated to be 7cg = l-5±o'?- Note that the 
estimated susceptibility exponent is significantly smaller 
than that of the standard 3D Ising EA model, 7 ~ 4 [Id]. 

If we combine the present estimate of 7cg with the 
estimate of Pcg from the off-equilibrium simulation of 
Ref. § and use the scaling relations, various chiral-glass 
exponents can be estimated to be a ~ —1.7, /3cg — 1-1, 
7cg — 1-5, vcg — 1-2 and t] C g — 0.8. The dy- 
namical exponent is estimated to be zgg — 4.7 by us- 
ing the estimated value of A and the scaling relation 
A = /Jcg/^cg^cg- While the dynamical exponent zgg 
comes rather close to the z of the 3D EA model, the 
obtained static exponents differ significantly from those 
of the 3D Ising EA model ~ 0.55, 7 ~ 4.0, v ~ 1.7 
and r\ ~ —0.35 [Id], suggesting that the universality 
class of the chiral-glass transition of the 3D Heisenberg 
spin glass differs from that of the standard 3D Ising spin 
glass. Possible long-range and/or many-body nature of 
the chirality-chirality interaction might be the cause of 
this deviation. 
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FIG. 2. Log-log plots of the time dependence of the equi- 
librium spin (a) and chirality (b) autocorrelation functions 
at several temperatures. The lattice size is L = 16 averaged 
over 64 samples. In (b), the best straight-line fit is obtained 
at T/J = 0.16, represented by the broken line. 
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FIG. 3. Temperature and size dependence of the G param- 
eter of the chirality defined by Eq. (9). The inset represents 
a log-log plot of the reduced chiral-glass susceptibility versus 
the reduced temperature t =| (T — Tcg)/Tcg |- 

Further evidence of a phase transition is obtained from 
the behavior of the G parameter of the chirality defined 
by 



G(L) 



H) 2 } H)? 
M)] - M)\ 2 



(9) 



While this quantity was originally introduced to repre- 
sent the non-self-averaging character of the system [pTj , 
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Bokil et al argued that it was not necessarily so p2[ . 
Still, a crossing of G(L) for different L, if it occurs, can 
be used to identify the transition 12 1. As shown in Fig. 3, 
for T > Tcg G(L) decreases with increasing L tending 
to zero, while for T < Tcg it tends to increase with L, 
thus lending further support to the existence of a phase 
transition at T = Tcg- 
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FIG. 4. Chiral-overlap distribution function below Tcg- 
The temperature is T/J — 0.1. 



We note that, if P{q x ) has a nontrivial structure as sug- 
gested from Fig. 4, the denominator of Eq. (9) should 
remain nonzero at < T < Tcg- Then, our data of G(L) 
in Fig. 3 indicates that the chiral-glass state is non-self- 
averaging. 

In summary, spin-glass and chiral-glass orderings of 
the 3D Heisenberg SG are studied by Monte Carlo sim- 
ulations. Our observation both on statics and dynamics 
strongly suggests the existence of a stable chiral-glass 
phase at finite temperatures without the conventional 
spin-glass order. This fact strengthens the plausibility 
of the chirality mechanism for experimentally observed 
spin- glass transitions. The nature of the chiral-glass or- 
dered state as well as of the critical phenomena are dif- 
ferent from those of the 3D Ising SG, with strong similar- 
ities to the system showing the one-step- like RSB, while 
its exact nature and physical origin have remained to be 
understood. 

The numerical calculation was performed on the Fu- 
jitsu VPP500 at the supercomputer center, ISSP, Univer- 
sity of Tokyo, on the HITACHI SR-2201 at the supercom- 
puter center, University of Tokyo, and on the CP-PACS 
computer at the Center for Computational Physics, Uni- 
versity of Tsukuba. 



In Fig. 4, we display the distribution function of the 
chiral-overlap defined by P (q') = [{S(q x — ql))] calcu- 
lated at a temperature T/J — 0.1, well below the chiral- 
glass transition temperature. The shape of the calculated 
P{q x ) is somewhat different from the one observed in the 
standard Ising-like models such as the 3D EA model or 
the mean- field SK model. As usual, P(q x ) has standard 
"side-peaks" corresponding to the Edwards- Anderson or- 
der parameter ±Qqq, which grow and sharpen with in- 
creasing L. The extracted value of ±<7cg coincides with 
that evaluated from the relaxation of C x {t). In addi- 
tion to the side peaks, an unexpected "central peak" 
at q x = shows up for larger T, which also grows and 
sharpens with increasing L. This latter aspect, i.e., the 
existence of a central peak, is a peculiar feature of the 
chiral-glass ordered state never observed in the EA or 
SK models. Since we do not find any sign of a first-order 
transition such as a discontinuity in the energy, the spe- 
cific heat nor the order parameter q^Q, this feature is 
likely to be related to a nontrivial structure in the phase 
space associated with the chirality. We note that this 
peculiar feature is reminiscent of the behavior charac- 
teristic of some mean-field models showing the so-called 
one-step RSB [lb]. Indeed, the existence of a negative 
dip in the Binder ratio gcG & n d the absence of a cross- 
ing of <7cg at <7CG > are consistent with the occur- 
rence of such one-step-like RSB Q. Our data of P(q x ) 
are also not incompatible with the existence of a con- 
tinuous plateau between [— 9cG'?cg] m addition to the 
delta-function peaks. According to the chirality mech- 
anism, such novel one-step-like RSB should be realized 
in the spin ordering of real Hciscnberg-like spin glasses. 
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